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We investigate the dimensional crossover of scaling properties of avalanches (domain-wall jumps) 
in a single-interface model, used for the description of Barkhausen noise in disordered magnets. 
By varying the transverse aspect ratio A = L y /L x of simulated samples, the system dimensionality 
changes from two to three. We find that perturbing away from d = 2 is a relevant field. The exponent 
r characterizing the power-law scaling of avalanche distributions varies between 1.06(1) for d — 2 
and 1.275(15) for d = 3, according to a crossover function f{x), x = {L'^/A, with <j> = 0.95(3). 
We discuss the possible relevance of our results to the interpretation of thin-film measurements 
of Barkhausen noise. We also study the probability distributions of interface roughness, sampled 
among successive equilibrium configurations in the Barkhausen noise regime. Attempts to fit our 
data to the class of universality distributions associated to l// a noise give a ~ 1 — 1.1 for d = 2 
and 3 (provided that suitable boundary conditions are used in the latter case). 

PACS numbers: 05.40.-a, 75.60.Ej, 05.65.+b, 75.50.Lk, 75.40.Mg 
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I. INTRODUCTION 

The Barkhausen effect [l| has long been known in mag- 
netism, and reflects the dynamics of domain- wall motion 
in the central part of the hysteresis cycle in ferromagnetic 
materials. The intermittent character which is a cen- 
tral feature of Barkhausen "noise" (BN) comes to light 
already in the original experimental setup. By wrap- 
ping a coil around a sample and ramping an external 
magnetic field at a suitable driving rate, one can detect 
well-separated voltage pulses across the coil, which are 
induced by sudden changes of magnetic flux. These in 
turn result from the microscopic realignment of groups 
of magnetic moments parallel to the field, i.e. domain- 
wall motion. For slow driving rates, the integral of the 
voltage amplitude of a given pulse over time is propor- 
tional to the change in sample magnetization, thus giv- 
ing a measure of the number of spins overturned in that 
particular event, or "avalanche size", to recall the termi- 
nology frequently used in the study of intermittent phe- 
nomena. Modern experimental techniques allow direct 
observation of domain-wall motion via magneto-optical 
Kerr effect measurements 0, Q , which demands use of a 
thin-film sample geometry. 

Early proposals for theoretical modelling of BN are re- 
viewed in Refs.0,0, whose authors formulate a Langevin 
description via Fokker-Planck equations. More recently, 
theoretical interest in the description of the statistical 
properties of BN has been rekindled, as attempts have 
been made to establish connections with general theories 
of non-equilibrium phase transitions and noi se p henom- 
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. Here, we shall be concerned with two such connec- 
tions. The first, motivated by the thin-film results just 
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alluded to, is the crossover between universality classes 
[these latter to be properly defined in the context] as one 
varies the spatial dimensionality of samples; secondly, we 
extend recently-developed concepts of universality of dis- 
tribution functions for l/f a noise M M M HH to the 
scaling properties of domain-wall roughness in BN. 

Experimentally, double-logarithmic plots of frequency 
of avalanche occurrence, P(s), against size s turn out to 
produce unequivocally straight sections, P(s) ~ s~ T , of- 
ten spanning 3-4 orders of magnitude, before dropping 
to zero for larger sizes 0, 01- Such power-law distri- 
bution of events has been associated to the concepts of 
self-organized criticality [6, 7, 25, 26], although other re- 
searchers argue that this in fact reflects proximity to a 
standard second-order critical point, together with an un- 
usually broad critical region in parameter space pi. ll2T|. 
Whatever the interpretation, a power-law decay shows 
up in assorted models used for numerical simulation of 
BN, both those based on the motion of a single inter- 
face in a disordered medium 0, 0, E3, El and those which 
adopt a picture of nucleation of multiple domains in a 
random-field Ising system 

Analogy with the well-established scaling theory of 
equilibrium phase transitions suggests that, in this case 
of a non-equilibrium phenomenon, the search for distinct 
universality classes may lead to a better understanding 
of the basic mechanisms involved. In Ref . HH, experimen- 
tal measurements of the exponent r of avalanche dis- 
tributions for several soft ferromagnetic materials were 
found to separate into two distinct groups, namely t = 
1.50 ± 0.05 (polycrystalline Fe-Si and partially cristal- 
ized amorphous alloys) and r = 1.27 ±0.03 (amorphous 
alloys under stress). It was then proposed that BN for 
each group of materials listed above belongs to a dif- 
ferent universality class of non-equilibrium phase transi- 
tions. While the value of t is a fairly plausible indicator 
of universality, or lack thereof, between different systems, 
many questions (prompted again by analogy with static 
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critical phenomena) still remain, such as to how many 
independent exponents there are, and, of particular in- 
terest here, what is the effect of s pac e dimensionality. 

Although sample shapes in Ref.Hslwere ribbon- like (30 
cm x 0.5 cm x 60 fim, to quote typical dimensions), this 
is far beyond the thin-film regime, for which thicknesses 
are of the order 5 — 100 nm (2, L3| Thus the behavior 
reported in Ref. is expected to be characteristic of 
fully three-dimensional objects. However, when consid- 
ering ever thinner samples as in Refs. 0, 0- dimensional 
crossover effects cannot be ruled out from the outset. 

In this work we use a single-interface model, originally 
introduced in Ref. yj for the description of BN. We re- 
call (though a detailed discussion will be deferred to Sec- 
tion that, for a fixed space dimensionality d = 2 or 
3, numerical values of e.g. the exponent r have been 
found to differ between nucleation and single-interface 
models (or even between distinct formulations of the lat- 
ter). Here we work under the assumption that the gen- 
eral features of dimensional crossover to be uncovered are 
model- independent, similarly to the presence of power- 
law avalanche distributions. The same assumption is ex- 
pected to hold as regards the roughness distributions to 
be investigated in Section Hvl 



II. MODEL AND CALCULATIONAL METHOD 

Here we shall use the single-interface model introduced 
in Ref. Q. We restrict ourselves to the adiabatic limit of 
a very slow driving rate, meaning that avalanches are 
regarded as instantaneous (occurring at a fixed value of 
the external field). Many experimental setups can be 
properly described in this approximation |^l7l llfl[lfl . ll£j . 

Simulations are performed on an L x xL y x oo geometry, 
with the interface motion set along the infinite direction. 
The interface at time t is described by its height hi = 
h{x,y,t), where (x,y) is the projection of site i over the 
cross-section. No overhangs are allowed, so h(x,y,t) is 
single- valued. Each element i of the interface experiences 
a force of the form: 



fi = u(x, y, hi) + k 



hij (i) - hi 



He , (1) 



where 



tension) type. In this term, £j(i) is the position of the j- 
th nearest neighbor of site i. The tendency of this term 
is to minimize height differences among interface sites: 
higher (lower) interface elements experience a negative 
(positive) force from their neighboring elements. The 
force constant k gives the intensity of the elastic cou- 
pling, and is taken here as the unit for /. We assume 
the boundary conditions to be periodic along x and free 
along y, so sites at y = and y — L y represent the film's 
free surfaces and have only three neighbors on the xy 
plane (except in the monolayer case L y = 1 which is the 
two-dimensional limit, where all interface sites have two 
neighbors). The last term is the effective driving force, 
resulting from the applied uniform external field H and a 
demagnetizing field which is taken to be proportional to 
M = (l/L x L y ) Ylt=i v hi, the magnetization (per site) of 
the previously flipped spins for a lattice of transverse area 
L x L y . For actual magnetic samples, the demagnetizing 
field is not necessarily uniform along the sample, as im- 
plied in the above expression; even when it is (e.g. for a 
uniformly magnetized ellipsoid), rj would depend on the 
system's aspect ratio. Therefore, our approach amounts 
to a simplification, which is nevertheless expected to ca 
ture the essential aspects of the problem. See Ref. 
for a detailed discussion. Here we use R — 5.0, k = 1, 
r\ = 0.05, values for which fairly broad distributions of 
avalanche sizes and roughness are obtained. 

We start the simulation with a flat wall. All spins 
above it are unflipped. The applied field H is set close 
to the saturation value of the effective field H e , in order 
to minimize transient effects. The force /, is then cal- 
culated for each unflipped site along the interface, and 
each spin at a site with fi > flips, causing the inter- 
face to move up one step. The magnetization is updated, 
and this process continues, with as many sweeps of the 
whole lattice as necessary, until fi < for all sites, when 
the interface comes to a halt. The external field is then 
increased by the minimum amount needed to bring the 
most weakly pinned element to motion. The avalanche 
size corresponds to the number of spins flipped between 
two interface stops. 



III. SCALING OF AVALANCHE 
DISTRIBUTIONS AND DIMENSIONAL 
CROSSOVER 



H e = H - T)M . (2) 

The first term on the RHS of represents the pinning 
force, u, and brings quenched disorder into the model by 
being chosen randomly, for each lattice site fl = (x, y, hi), 
from a Gaussian distribution of zero mean and standard 
deviation R. Large negative values of u lead to local 
elements where the interface will tend to be pinned, as 
described in the simulation procedure below. The second 
term corresponds to a cooperative interaction among in- 
terface elements, assumed here to be of elastic (surface 



We have collected avalanche histograms for varying L x , 
L y , in a such a way that the number of interface sites 
L x L y varies between 800 and 80,000. The aspect ratio 
A = L y /L x was varied between essentially zero (d = 2, 
one-dimensional interface) and unity (d = 3, square in- 
terface). For each L x , L y we generated 10 5 avalanches. 
Although it may take 10 2 — 10 3 avalanches for a steady- 
state regime to be reached (as measured by the stabi- 
lization of H e against external field H, apart from small 
fluctuations), we have checked that the only distortion 
introduced by the transient on avalanche size statistics 
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Figure 1: Double- logarithmic plot of avalanche size distribu- 
tion for L x = 60, L y = 30 (circles). Full line is a fit to the form 
Eq. JIJ, for which the optimal parameters are: r = 1.226(6), 
5 = 3.6(1), sO = 5.97(5) x 10 4 . 

is the one large event occurring at the very start, i.e. on 
departure from the initial, flat-interface, configuration. 
We cannot guarantee this to be so when roughness is 
the quantity under investigation, thus data in Section Hvl 
have been collected only under steady-state conditions. 

The probability distribution P(s) for avalanche size s 
is expected to behave as 

= .,-/(£) , (3) 

where so is a cutoff related (in experiment) to domain 
size and/or demagnetization effect s HTl ll^|. and (in sim- 
ulations) to finite-lattice effects [?lll0l|. or proximity to a 
critical point @. The specific shape of the function f(x) 
has been debated. While a simple exponential has often 
been used, either phenomenologically or fin 

some special cases) backed by theoretical arguments @, 
Gaussian fits, f(x) ~ exp(— s 2 / s o) nave been pro- 
posed as well. Going one step further, and at the same 
time trying to keep the number of fitted parameters to 
a minimum, here we shall follow Refs 113, and fit our 
data to a stretched exponential: 

P( s ) = A s- T e- {s/so)5 . (4) 

Apart from the overall normalization factor Aq, one then 
has three free quantities to fit, which has proved enough 
for our purposes. A typical example is displayed in Fig.Q] 
For the roughly sixty {L x , L y } sets investigated here, the 
fitted value of 6 usually falls in the interval 2.4 — 3.5, with 
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Figure 2: Effective exponent r from fits of simulation data 
to Eq. JIJ, against inverse cross-sectional interface area 
(LxLy) -1 . Bottom (full triangles linked by full line): mono- 
layer, L y = 1. Top (Full squares linked by full line): d — 3, 
L x = L y . Intermediate curves (open symbols connected by 
dashed lines): from bottom to top, aspect ratio A = 0.005, 
0.01, 0.02, 0.1 . 



half a dozen cases slightly above that. This is broadly in 
line with 6 = 2.32(6) quoted in Ref. Though one 
might argue that a universal value should hold for this 
exponent, we feel that our results are not accurate enough 
either to prove or disprove such hypothesis. 

On the other hand, the exponent r which is of central 
interest here displays systematic variations both against 
lattice dimensions and aspect ratio. These are displayed 
in Fig. [51 The overall picture strongly suggests a sys- 
tematic crossover towards three-dimensional behavior for 
any fixed (finite) aspect ratio. In addition to this, finite- 
lattice effects are present as well. 

In order to gain a quantitative understanding of this, 
we recall general ideas of finite-size HUHD] and crossover 
scaling |3ll l32l |. Crossover phenomena reflect the com- 
petition between different types of (pseudo-) critical be- 
havior in the vicinity of a multicritical point, at which 
several characteristic lengths diverge (i.e. their associ- 
ated fields approach their respective critical values). See 
e.g. Section VILA of Ref. £l] for an illustration of the 
well-known case of thermal-geometric crossover in dilute 
magnets near the percolation point. Closer analogy with 
the present case is found in the discussion of dimensional 
crossover in "layer" magnets jH HH (see fi gure 42 in 
Ref. Eul) . Finite-size scaling can be seen as a particu- 
lar instance of crossover, in which the system's inverse 
finite size L _1 is an additional relevant field, driving it 



4 



away from the true criticality which occurs only in the 
thermodynamic limit [2<1 Ojij, 03] . We now turn to the 
simulational results exhibited in Fig. |2 Upon increasing 
sample size, finite-size effects are reduced and the trends 
followed by the respective exponent estimates differ, de- 
pending on whether (i) the aspect ratio A is kept fixed, no 
matter how small its value (all curves except the lower 
one), or (ii) L y = 1 is fixed instead (monolayer, lower 
curve). The latter corresponds to ever-decreasing A as 
L x increases, and in the L x — » oo limit is expected to 
reflect true two-dimensional behavior. Considering, for 
instance, the curve for A = 0.005, for L x not very large 
one has exponent estimates closely resembling those of a 
monolayer with a similar cross-sectional "area". However, 
the trend shown implies that this is in fact an apparent 
behavior, which arises only as long as finite-size effects 
(represented by L' 1 ) are more important than the sys- 
tem's three-dimensional character (represented by a non- 
zero value of A). Increasing L x while keeping A constant, 
the relative importance of these quantities is eventually 
reversed, and one crosses over to three-dimensional be- 
havior. These remarks are translated into quantitative 
statements, as follows. Since the inverse finite size L~ x 
and the interface aspect ratio A are both relevant fields, 
which drive the system away from true two-dimensional 
behavior, a plausible ansatz for the crossover variable is 
x = (L^ 1 )^ I A. where tj> is a crossover exponent, to be 
determined pul Thus, for any fixed A / 0, and 

L x — » oo (x — » 0), three-dimensional features must dom- 
inate, while for x ^> 1 two-dimensional behavior (with 
finite-size corrections) will take over. 

Furthermore, the data of Fig. must collapse on the 
same curve, when plotted against x. This latter state- 
ment gives the operational procedure for determination 
of <f>. By recalling that the horizontal axis variable in 
Fig. can be written as (L x L y )- 1 = x 2 ^ A 2 ^~ l , and 
that constant-^ curves are further away from the verti- 
cal axis the larger A is, one sees that, in order for those 
curves to collapse one must have (p < 2; consideration of 
the numerical values involved shows that, in fact, (f> < 1 
is needed. The choice of <f> — 0.95(3), with the corre- 
sponding results exhibited in Fig. EJ reflects the range of 
<j> for which the collapse of all data for non-zero A (onto 
the left-hand side of the diagram) is visually deemed to 
be best. The fact that the d = 2 data are segregated to- 
wards the right is to be expected in this context, as they 
belong to a different side of the crossover (x > 1) where 
three-dimensional effects vanish. 

The intersection of the scaling curve with the vertical 
axis provides an estimate of the three-dimensional scal- 
ing exponent. The result is r(d — 3) = 1.27(1). An ad 
hoc parabolic fit against 1/L X , including only data for 
A = 1 and L x < 160, gives r(d = 3) = 1.28(1), where 
the largest contribution to the uncertainty comes from 
spread between extrapolations that either do or do not 
include small-lattice data (L x — 30, 40). As we have 
seven lattice sizes available for A = 1, use of parabolic 
fits already gives us at most four degrees of freedom. 
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Figure 3: Effective exponent r from fits of simulation data 
to Eq. (QJl, against crossover variable x = (L" 1 )^ /A, with 
4> — 0.95. Key to symbols is same as in caption to Fig. [2] . 
Full line is a fourth-degree fitting curve with an exponential 
tail near x = 0. Inset shows details of main figure close to the 
vertical axis. 



Given the limited range of data available, attempting to 
improve estimates by including higher-order corrections 
would therefore not seem justifiable. Since neither ex- 
trapolation procedure appears to be obviously superior 
to the other, our final quote encompasses both results: 
r(d = 3) = 1.275(15). 

In Ref.Q, r(d = 3) = 1.13(2) is quoted for L x = 40, 
while our corresponding result is 1.21(1). This difference 
arises mainly from distinct fitting procedures, in partic- 
ular those authors' apparent use of a fixed 8 = 1 for the 
cutoff. Indeed, by keeping 6 = 1, the best fit of our data 
is for r(d = 3) = 1.12(4), though at the cost of increasing 
the x P er degree of freedom by one order of magnitude 
compared with the variable-(5 fitting scheme. 

A similar parabolic fit of the two-dimensional data 
gathered on the right-hand corner of Fig. (using 4000 < 
L x < 16000), gives r(d = 2) = 1.06(1) . This is 
broadly in accord with the result of Ref. Q, namely 
r(d = 2) = 1.00(1) (using 500 < L x < 5000). 



IV. ROUGHNESS DISTRIBUTIONS AND l/f a 
NOISE 

In this section, we consider only two-dimensional 
(L x = 1) and fully three-dimensional (L x = L y ) sys- 
tems. We have collected data on interface roughness, in 
order to analyze them in the context of universal flue- 
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tuation distributions. An important reason for interest 
in universal distributions is that they have no adjustable 
parameters The fact that a given prop- 

erty of a system behaves according to one of such dis- 
tributions is then expected to indicate the universality 
class to which it belongs. For non-equilibrium problems 
such as is the case here, the task of connecting exponents 
(or features of distributions) to universality classes is far 
from accomplished (see e.g. the discussion in Ref . I23T1 . 
Accordingly, our purpose here is simply to identify what 
universal distribution (if any) is followed by the scaled 
interface roughness for the current model. 

In order to make contact with previous work on inter- 
face roughness at depinning transitions, we recall some 
basic ideas. 

At the end of each avalanche, we measured the rough- 
ness 1Z of the instantaneous interface configuration at 
time t in the usual way, as the rms fluctuation of the 
interface height: 



n(t) = 



L X L V 



(LxLy)- 1 £ (hi(t)-h(t))' 



1/2 



(5) 



where h(t) is the average interface height at t. As the 
avalanches progress, one gets a sampling of successive 
equilibrium configurations; the ensemble of such config- 
urations yields a distribution of the relative frequency 
of occurrence of 1Z. In order to get clean distributions, 
we have seen that the number of events considered must 
be O(10 6 ), i.e. one order of magnitude larger than the 
samples used for the analysis of size distributions in Sec- 
tion IfflJ We have used only steady-state data. i.e. after 
the stabilization of H e of Eq. (J2J against external field H . 
This way, our implicit assumption that successive inter- 
face configurations are stochastically independent gains 
plausibility. Similar ideas were invoked in Ref. to 
justify a mapping of the steady state of a deposition- 
evaporation model onto a random-walk problem. 

The roughness exponent C is related to the finite-size 
scaling of the first moment of the distribution 



(n L ) ~ l< 



(6) 



where the angular brackets stand for averages over the 
ensemble of successive interface configurations of an in- 
terface with transverse dimension L. We have estimated 
C from power-law fits to our data. For d = 2 (one- 
dimensional interface) using 400 < L x < 1200 we get 
((2d) = 1.24(1), which compares well with the usually 
accepted £ ~ 1.25 for the quenched Edwards- Wilkinson 
universality class El Hi H3- For d = 3 (two- 
dimensional interface) we used 30 < L x = L y < 80 
and two alternative sets of boundary conditions, namely 
mixed (MBC) i.e. free along x and periodic along y, as 
described in Section |nl and periodic along both x and y 
(PBC). The results are, respectively, ((d = 3, MBC) = 
0.87(1) and ((d = 3, PBC) = 0.71(1). While the lat- 
ter value is not far from C — 0.75 for the corresponding 



quenched Edwards- Wilkinson model [3(1 l39| . the former 
seems difficult to relate to existing results. 

At this point one might invoke universality ideas and 
claim that the difference between three-dimensional es- 
timates must be a finite-size effect. Investigating this 
directly, e.g. by performing simulations for larger L, 
would be straightforward but time-consuming. However, 
as we shall see in what follows, this issue can be ad- 
dressed more efficiently by returning to our main goal 
in this Section, i.e. by looking at the full distributions, 
rather than examining only selected moments (as is the 
case of the scaling for extraction of Q . This is because 
mounting evidence indicates that, in general, width dis- 
tributions decouple into the product of a single size- 
dependent scale and a universal (size-independent) scal- 
ing function IHIH that is 

P{K) = $(K/a) , a 2 = ((U (K)) 2 ) . (7) 

Therefore, if the functional form $(z) varies depending 
e.g. on whether mixed or periodic boundary conditions 
are used, one is on safer grounds to assume that this 
reflects differences in the respective universality class. 
In this context, it must be recalled that dependence of 
scaling quantities on boundary conditions is an often- 
encountered feature when dealing with fluctuation phe- 
nomena m,mi- One illustration of this, which is of great 
relevance here, is that one of the requirements for a class 
of universal (time) fluctuation distributions to hold is 
that they be periodic in time [12, 0,0] ■ 

In our case, this feature is replaced by periodic bound- 
ary conditions in space (the same reasonin g w as used for 
the deposition-evaporation model of Ref. I34R. Strictly 
speaking, this time-space correspondence is only true 
for the one-dimensional interfaces of the two-dimensional 
version of our model. However, the mixed boundary con- 
ditions (MBC) used in Sec. IIIII suggest that an extension 
of our investigation to the three-dimensional case may 
not be unjustified. As regards fully periodic boundary 
conditions (PBC), even though it is not obvious that the 
analogy can be pushed that far, we decided to analyse 
the respective data for completeness. 

We analysed the roughness distributions for d = 2 and 
for the fully three-dimensional case (i.e. aspect ratio 
L y /L x = 1), the latter both with MBC and PBC. As 
shown below, finite transverse dimensions are of negligi- 
ble import as far as the scaling functions <&(z) are con- 
cerned, thus giving further support to the assumption 
that Eq. © holds. 

We have compared our results against the family of 
roughness distributions for periodic 1//" noise, described 
in Ref.m- The roughness of a time signal, and its connec- 
tions via Fourier transform with the frequency spectrum 
(and thus with the respective a), are described at length 
in Section II of that reference. The idea of fitting spa- 
tial roughness data to the roughness of time signals (e.g. 
random- walks) is well-known in the study of interface 
fluctuations ■ The new feature here is that we have at 
our disposal a family of distributions whose shape varies 
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Figure 4: Scaled probability distributions $(y) in d — 2 for 
y defined in Eq. JjjJ. Crosses: = 400; circles: L x = 1000. 
Lines are roughness distributions for l/f a noise given in 
Ref.El , with a = 1.15 (full) and 1 (dashed). 



Figure 5: Scaled probability distributions $(y) in d — 3 with 
MBC, for ?/ defined in Eq. ©. Crosses: L = 30; circles: 
L = 80. Lines are roughness distributions for l/f a noise 
given in Ref. , with a = 1.15 (full) and 1 (dashed). 



smoothly against the single basic parameter a (but, apart 
from that, are strictly parameter-independent). We ex- 
pect that the value of a which best fits our data should 
be connected to the underlying universality class of thew 
avalanche model used here. We recall that ours is a 
phenomenological study, since at present very little is 
known regarding the (physical) causal relationships be- 
tween fluctuation distributions and universality classes 
for out-of-equilibrium transitions (53, EI E3, ■ 

The quantity to be compared against a universal form 
is the distribution of the deviation from the average, 
scaled by the standard deviation (this is termed scaling 
by the variance in Ref. 1 231) : 



Our data are depicted respectively in Fig.^lfor d = 2, 
Fig. for d = 3 with MBC, and Fig. El for d = 3 with 
PBC. In each figure, we have plotted data corresponding 
to only two values of L ("small" and "large"), in order 
to illustrate that the i-dependence of the scaled distri- 
bution is negligible, without cluttering the diagram with 
intermediate-L results. In all figures, the analytic curve 
corresponding to a = 1 is shown as a dashed line, in order 
to ease comparison between different cases. In general, 
the scatter of simulational data at the low-roughness end 
of the distribution is much smaller than that at the high 
end, thus we might give more weight to the former re- 
gion when judging the quality of fit. Using this criterion, 
one sees that the d = 2 data appear to be closer to the 



a = 1.15 curve than to that for a = 1, while the situation 
is reversed in the d = 3 case with MBC. Furthermore, for 
d = 3 with PBC the best fit is undoubtedly for a < 1. 
For this latter case, we also display a Gaussian distribu- 
tion, which corresponds to a — 0.5 j23|. We have tried to 
quantify the above remarks, by investigating the behav- 
ior of the x 2 P er degree of freedom (x 2 /d.o.f.) for fits of 
our data to the analytical distributions, against varying 
a. We have included only data for which $(y) > 10~ 3 , 
i.e. those displayed in Figs. 31 El an d El Though, in 
principle, the distributions can be evaluated in closed 
form [23], we ran into serious numerical problems for 
a < 0.9 in the region y < — 1. Fortunately, as shown 
in Fig. E| this does not matter much as long as the d = 2 
and d = 3 (MBC) cases are concerned, because the re- 
spective x 2 /d.o.f. clearly exhibit minima located slightly 
above a = 1. On the other hand, this means that for 
d = 3 with PBC we were not able to follow the trend 
shown in the Figure into the region a < 0.9, where it is 
clear that a minimum of the corresponding x 2 /d.o.f. must 
be located. On the other hand, the Gaussian (a = 0.5) 
distribution shown in Fig. El evidently overshoots the de- 
sired corrections, so one can be sure that the best fit will 
be in the interval 0.5 < a < 0.9. 



V. DISCUSSION AND CONCLUSIONS 

The data displayed in Section ITTT| show that, in the in- 
terface model introduced in Ref. Q, the power-law behav- 
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Figure 6: Scaled probability distributions <3?(y) in d — 3 with 
PBC for y defined in Eq. ©. Crosses: L — 30; circles: L — 
80. Lines are roughness distributions for l/f°" noise given 
in Ref. 23; , with a = 0.9 (full), 1 (dashed) and 0.5 (short- 
dashed). 
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Figure 7: x 2 P er degree of freedom (x 2 /d.o.f.) for fits of sim- 
ulation data to analytical forms of 1/ f a distributions, against 
a. Triangles: d = 2, L x = 400; Squares: d = 3 MBC, L = 40; 
hexagons, d = 3 PBC, L = 40. 



ior of avalanche statistics is characterized by an effective 
exponent, which varies continuously both with the inter- 
face's linear dimensions and aspect ratio. By means of a 
finite-size and crossover analysis, we have demonstrated 
that such continuous variation in fact reflects crossover 
towards three-dimensional behavior, for any non-zero as- 
pect ratio. 

The implications of this for the interpretation of ex- 
perimental BN results in thin films must be worked out 
carefully. In Ref.0, the value r = 1.1 is given for Fe films. 
The authors of Ref.0quote r ~ 1.33 for Co films and con- 
clude that their setup is a two-dimensional realization of 
the (single-interface) model of Ref. for which r = 4/3 
in d = 2 and 3/2 in d = 3 [2 El 0n the other hand, 
as seen above, the single-interface model considered here 
gives t = 1.06(1) in d = 2 and 1.275(15) in d = 3. Thus 
one might interpret both sets of experimental results as 
reflecting a crossover towards three-dimensional behav- 
ior. As regards the specific features of the experimental 
investigations, one must ask: (i) which, if any, of the 
two models applies to the corresponding microscopic de- 
scription, and (ii) how far along, quantitatively, is the 
dimensional crossover for the thin-film geometries used. 

While the visual evidence displayed in Ref. is con- 
vincing proof that a single-interface picture applies in 
that case, for a definite answer to (i) one must look at the 
differences between the models. As far as the power-law 
scaling of avalanche distributions is concerned, the model 
of Ref. differs from the one considered here by the intro- 
duction of a non-local kernel due to dipolar interactions. 
The value r = 4/3 quoted in Ref. relies on assuming 
that the form taken by this kernel in momentum space 
is ~ <7 M , (q = wavevector) with fi = 1 [i^J. In d = 3 
the same theory gives r = 3/2 for // = 1, and r = 5/4 
for /i = 2 [Igj. Both values have been found to good 
approximation in experiments on fully three-dimensional 
systems (thus defining distinct universality classes), as 
recalled in the Introduction 0|. For the thin-film cases 
it is not clear, without a detailed analysis of the spe- 
cific materials involved, whether the non-local kernel is 
of sufficient import to drive avalanche scaling towards the 
d = 2 behavior predicted for the model of Ref. 0. 

Turning to question (ii) , recall that the evolution of an 
interface along a 400 x 320 /zm 2 area of a 25-nm film is 
shown in Ref. Translating to the language of Section 
IIII1 this would correspond to a transverse aspect ratio 
A = 25nm/320/im ~8x 10~ 5 (this in an upper bound, 
as the film's transverse dimensions are likely to be larger 
than the area shown). Though we do not think that 
the model results depicted in Fig. are quantitatively 
accurate enough, one must keep in mind the possibil- 
ity that the effective experimental behavior still is very 
close to the two-dimensional limit. Indeed, the simula- 
tional curve for ^4 = 5xl0~ 3 already shows a value of r 
rather close to that for the two-dimensional case, along 
an extended portion of the Figure. The surest way to 
settle this matter would be by performing a series of ex- 
periments on films of the same composition and varying 
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thicknesses, in order to produce a full picture of the di- 
mensional crossover. We hope experimentalists will be 
motivated by the present results. 

As regards the search for universal roughness distri- 
butions in Section IWl for now we quote (from Figured 
and the associated remarks) a ~ 1.05 (d = 2); a ~ 1.15 
(d = 3, MBC) and 0.5 < a < 0.9 (d = 3, PBC). It 
thus appears that the boundary conditions do have sig- 
nificant influence in this context, a fact which remains to 
be more fully understood. Although the search for the 
physical origins of 1// noise 0] is clearly of great inter- 
est, it appears that, at least for non-equilibrium phenom- 
ena as is the case here, we are still at a very preliminary 
stage. Again, it is hoped that the present results will mo- 
tivate further research. Measurements of the roughness 
distributions for alternative models of BN p, Ll| would 



be a natural extension of the this work, in order to check 
whether the above-quoted values of a are indeed univer- 
sal within this subset of avalanche models. 
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